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Abstract 

Deformation quantization (sometimes called phase-space quantization) 
is a formulation of quantum mechanics that is not usually taught to un- 
dergraduates. It is formally quite similar to classical mechanics: ordinary 
functions on phase space take the place of operators, but the functions 
are multiplied in an exotic way, using the ^-product. Here we attempt a 
brief, pedagogical discussion of deformation quantization, that is suitable 
for inclusion in an undergraduate course. 

1 Introduction 

Another way of doing quantum mechanics grew from pioneering works of Weyl, 
Wigner, Groenewold, Moyal, Baker, and others. For reviews, see pQ. Its im- 
portance as an autonomous formulation of quantum mechanics was appreciated 
in the 1970s; it was also understood that it could be viewed as a deformation 
of classical mechanics [2]. The formulation was therefore dubbed deformation 
quantization. 

The pedagogical review article [3] advocated that deformation quantization 



be included in graduate studies, and it is treated in some graduate textbooks 
(see for example). We believe undergraduate students would also benefit by 
learning something about deformation quantization, i.e. from an exposure to it. 
For that reason, we attempt here a brief, pedagogical discussion of deformation 
quantization, with upper-level undergraduates and their instructors in mind. 

To make our presentation as pedagogical as possible, we will restrict to the 
case of one degree of freedom, i.e., to the case of a single particle moving on the 
x-axis. We'll only treat pure quantum states, since the generalization to mixed 
states is straightforward. Lastly, although the Weyl ordering and correspond- 
ing Groenewold-Moyal ★-product determine just one of many different ways to 
do deformation quantization, we will only discuss the Weyl-Groenewold-Moyal 
case. 

It is important to point out, however, that deformation quantization is not 
just of pedagogical interest. Recently, it has been an active research topic 
in both physics and mathematics. Physicists studying string theory use the 
methods of deformation quantization because in certain conditions, strings live 
in spaces whose coordinates do not commute [5] |H] , much as the two quantum 
coordinates X, P of phase-space do not, according to ©. The mathematician 
Kontsevich's work on deformation quantization in 7j was part of the reason he 
was awarded the Fields medal, math's highest honor. 

2 Overview 

Consider a single particle, with position x and momentum p. Phase space is 
the two-dimensional space with coordinates (x,p). Each point of phase space 
specifies a classical state of the system, and as a state evolves in time, a point 
traces out a path in phase space. For example, a simple harmonic oscillator 
follows an elliptical trajectory, centered on the origin (x,p) = (0,0). 

If our knowledge of a system were imprecise, the state might be given as a 
probability distribution on phase space, perhaps a bump with its center at the 
most probable values (x,p). Classical dynamics could be done by following the 
evolution of this distribution. 

Can quantum mechanics be done in a similar way? The answer is yes, and 
the way is called deformation quantization 

By doing quantum mechanics this way, the introduction of abstract quantum 
states, and their Hilbert spaces, can be avoided. States are described instead 
by functions on phase space, as in classical mechanics. As a consequence, the 
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relation between quantum and classical mechanics may be understood better. 

Before sketching how deformation quantization is done, we need to empha- 
size that quantization of a classical system is not a unique procedure, no matter 
what formulation of quantum mechanics is used. In the operator approach, x 
and p are replaced by the corresponding operators, denoted X and P. Sup- 
pose we need to work with something like x 2 p in quantum mechanics, should 
we consider X 2 P1 Or XPX, for example? The ambiguity can be reduced 
by demanding that the operators constructed be Hermitian, but that does not 
eliminate the choice completely. 

This operator-ordering ambiguity is not a new problem, special to deforma- 
tion quantization, since it is part of the usual operator approach to quantum 
mechanics. A choice must be made, so let's use the so-called Weyl ordering, 
and write 

ti{x 2 p) = - (X 2 P + XPX + PX 2 ) . (1) 
o 

Weyl ordering can be extended to functions on phase space by specifying how it 
works on all monomials x m p n , and then applying it to the Taylor expansions of 
functions. ■d(x m p n ) is simply the average of all possible orderings of m factors 
of X and n factors of P. 

i? is known as the Weyl map, taking functions on phase space to operators. 
We now have an operator $(/) that corresponds to a function / = f(x,p) on 
phase space. We'll call such an operator a Weyl operator. Multiplying two Weyl 
operators gives another one, so their algebra closes, as it must. 

Remarkably, one can prove a stronger statement. It is this result that makes 
deformation quantization possible. Groenewold showed that 

*(/)<%) = #(f*9) ■ (2) 

That is, multiplying two Weyl operators is equivalent to *- multiplying the cor- 
responding phase-space functions, and then applying The ^-product (pro- 
nounced star-product) takes the form 

f{x,p)*g(x,p) = f(x,p) expjy ( dxdp - dpdx )} g{x,p) . (3) 

Here d x = etc., and the arrows indicate the directions in which the deriva- 
tives act. The exponential is to be understood using the series expansion 



Since it is directly related to the product of Weyl operators, the ^-product is 
non-commutative and associative, as the product of operators is. It is a strange- 
looking product precisely because it must mimic the product of operators. 

Eqn. (01 is important because it suggests that one might be able to avoid 
constructing operators from functions on phase space and just work with the 
functions directly, as long as they are multiplied using the ^-product. This 
is exactly what is done in deformation quantization. Operator products are 
changed to ^-products, and the Weyl map $ is factored off, roughly speaking. 

To give the flavor of how it goes, it is easy to show from J3J that 

X-kp = xp + — (5) 

and 

p-kx = xp - — . (6) 

Then the *- commutator, or Moyal bracket, of x and p is 

= x -k p — p-kx — ih . (7) 
This result is consistent with the crucial canonical commutation relation 

[X, P] = ih. (8) 

To do quantum mechanics we need further ingredients, beyond the ^-product. 
The quantum state must be described. In the operator formulation, a pure quan- 
tum state is describable by a state vector \tp) . The most general type of quantum 
state is a mixed state, however, and it incorporates classical probabilities for dif- 
ferent pure states. It is a mixed quantum state that corresponds to the classical 
distribution on phase space mentioned in the second paragraph. To specify such 
a mixed state, the density matrix (sometimes called the state operator) p must 
be used (see [S] for a nice discussion). 

We'll nevertheless restrict this discussion to the case of a pure state if>, since 
it makes the presentation easier to follow. Then 

P = iv) m . (9) 

Generalization to mixed states is straightforward. 

Now we have an operator, and to do deformation quantization, we need a 
function on phase space that corresponds to it in the manner of J2J). That is, 
we need an operation W acting on operators and giving phase-space functions, 
that satisfies 

- / , (10) 
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for any / = f(x,p); we need W = $ . What works is 

W(G) = h Jdye-^ix + ^Glx-^) . (11) 

As with X, P vs. x,p, we will denote operators with upper-case symbols, to dis- 
tinguish them from phase-space functions. W(G) is called the Weyl transform 
(and sometimes the Weyl symbol) of the operator G. 

The Weyl transform of the density matrix p = \ip) (tp\, 

W(p) = h Jdye-^rix- ^)V(*+ y) , (12) 

is the central object in deformation quantization. After normalization, it is 
known as the Wigner function: 

>-^- 

It describes the quantum state of the system, and all observable probabilities 
can be calculated from it. 

This hints at a punchline: deformation quantization is the Weyl transform 
of quantum mechanics done with the density matrix, or state operator. 



3 Weyl Transform and Groenewold-Moyal Star 
Product 

In this section, some detail and proofs omitted in the previous section will be 
provided. It can be skipped in a first reading. 

First, consider the Weyl map of functions on phase space, like f(x,p). By 
expanding, it is not hard to convince oneself that 

ti((ax + bp) n ) = (aX + bP) n , (14) 

for all parameters a, b. It follows that 

d{e ax+bp ) = e aX+bp , (15) 

using Q. Now, the Taylor series of f(x,p) about (x,p) = (0,0) can be written 

as 

f(d a ,d b )e ax+b P\ a . b=0 = f(x,p) . (16) 
This combined with 115|) gives 

m = f(d a ,d b )e aX+bp \ a , b=0 , (17) 
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a useful general formula for Weyl operators. Slightly different versions of the 
same formula are 

= /M)e^ x+ ^ P k 6=0 = e^ x+ ^ p /(a,6)| a , fc=0 • (18) 

Another formula that one sees quite often is 

«?(/) - / drdadxdp f(x,p)e lT( - p - p ^+^ x -^ . (19) 

It follows from (|15|l and Fourier methods. The usual Fourier expression is 

f(x,p) = J drdo f{T,v)e i ° x+iTP , (20) 

where 



f(r,a) = ™ y dx'dp' f(x',p')e-^-*™ (21) 

is the Fourier transform of f(x,p). According to l|15l) . the Weyl map $ simply 
replaces e lax+lTp with e ^+^^. Making that replacement in Q20II21[I gives (|TT)|) . 

Equations 1)1711181111^1 are not so useful for calculating simple examples. They 
are, however, important for discussing the general properties of Weyl maps. For 
example, the -^--product is introduced because of the crucial property ©. In- 
order to prove it, l|18f) can be used. First, 

<?(/)<%) = f(a,b)eV« x ^ p e°* x ^ p g(a,b) L lM , 5=0 ■ (22) 
Now use the (simplified) Baker-Campbell-Hausdorff formula, 

e A e B = e ±[A,B] e A+B ^ (23) 

valid when the commutator [A, B] commutes with both A and B. Eqn. 
becomes 

mm = e^ +d ^ d ^ p f(a,b)ef $-°*-°>°*hM \ a , blSli=0 , (24) 



or 



= e^ x+d " p fia^e^^" 9 "-^ 9 ^ g(a,b) \ a , b=0 . (25) 
The result @ then follows. 

As mentioned above, the Groenewold-Moyal ^-product, defined in |J5J, is 
non-commutative, i.e. f * g ^ g * f ■ It is, however, associative 

(/i*/ 2 )*/ 3 = /i*(/ 2 */ 3 ) • (26) 
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It shares those properties with the product of operators, as demands. Simi- 
larly, it is easy to show that 

(a*/ 2 )* = / 2 **/r , (27) 

in agreement with the rule (-Fi-Fi)^ = f° r operators F\, F^. 

The exponent of the ^-product J2J indicates the most important property of 
deformation quantization: its intimate relation to classical physics. In classical 
mechanics, it is the Poisson bracket of functions on phase space, 

<*»>- t% ~ II = naA-w.), <*) 

that enters the dynamical equations [2] • In the operator formulation of quantum 
mechanics, it is the commutator [F, G] of operator observables F and G that is 
important. In deformation quantization, the ^-commutator 

[f,gl ee f*g-g*f (29) 

(recall eqn. 10) of functions / and g takes its place. The equation 

hm ![/,<?]* = {/,<?} (30) 
h^o iti 

encodes the relation between classical and quantum mechanics in deformation 
quantization. 

Now let us consider the Weyl transform W(G) of an operator G, defined 
by the property IjlOl) . and given explicitly by the formula (|llf) . Of course, here 
G = G(X, P) is an operator function of the operators X and P. 

First, let us argue for the formula I jllj l. From the development preceding 
H17|) . it is clear that we only need to show 

w( - e aX+hP) = e ax+bp _ ( 31 ) 

i.e., the inverse of i|15|) . This is not difficult, however, using the Baker-Campbell- 
Hausdorff formula (|23|l . and that e hp is a translation operator, i.e., e bp \x) — 
\x + ifib). 

Finally, we should note that there is an inverse analog of ©: 

W(G)*W(£) = W(GL) . (32) 

Because of fliup i . it works if the operators are Weyl operators, i.e., if G — $(g) 
and L = t9(£) for some phase-space functions g and I. They will be, however: 

X 2 P 2 X = i^{x){){x)d{p)d{p)d{x) = ${x*x*p*p*x) , (33) 

for example. 
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4 Wigner Function 



Now that the basics have been established, we can consider the description of 
quantum states and their evolution. As stated above, the central object is the 
Wigner function H12|) . 

In the operator method, one must determine the description of a state. Its 
state vector the description, is found by solving the Schrodinger equation 
4z\ip) = | ip) = \H\ip). In a similar way, the starting point in deformation 
quantization is the dynamical equation for the Wigner function p. From (|12|) . 
we find 



!-s/*'-~{<* + tI*»<*-t> 

+ {x+^L\i,){i,\x-^-)} . (34) 



Substituting the Schrodinger equation (and its adjoint) then gives 

% - 

= -j [H,p]i, ■ (35) 
in 

In the last step, we used (|32|) . and we have defined the Weyl transform H = 
W(H) of the operator Hamiltonian H . 

For stationary states, ^ = 0, so that 

[H, p]„ = . (36) 

Thus the Hamiltonian and Wigner function ★-commute. A stronger relation 
can be derived more directly. The Schrodinger equation simplifies for stationary 
states to H\ip) = E\ip), where E is the energy. This implies that Hp = Ep, 
which Weyl transforms to 

H*p = Ep . (37) 

In the next section, this simplified dynamical equation will allow us to solve for 
the Wigner function of the stationary states of the simple harmonic oscillator. 

Once the Wigner function is determined, how is it used? First of all, by 
(|12l) . the probability densities are 

V{x) = |?/>(x)| 2 = J dpp{x, P ) , 

V{p) = \^{P)\ 2 - / dx P (x lP ) . (38) 
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Clearly then, the Wigner function is normalized and real: J dxdpp — 1 and 

p = p- 

All observable expectation values can be calculated using the Wigner func- 
tion. The expectation value of an operator G is 

(G) = (V>|G|V> = Tr(pG) . (39) 

In deformation quantization, this translates into 

(G) = J dxdp p-kg , (40) 

where g = W(G). Roughly, one can think of the integral over phase space as 
the analog of the trace, and as discussed above, the star product takes the place 
of the operator product. The important cyclic property of a trace is encoded in 

J dx dp / * g = J dx dp g * / = Jdxdpfg. (41) 



5 Example: Simple Harmonic Oscillator 



The most common non-trivial example studied in physics, quantum or classical, 
is the simple harmonic oscillator (SHO). We will now treat the quantum SHO 
using deformation quantization, following 10 . Our results will illustrate some 
of the properties of the Wigner function in deformation quantization. 

Recall from above that H * p = Ep. Taking u = m = 1 for simplicity, the 
SHO Hamiltonian is H = \{P 2 + X 2 ). To use (EJ, we need to calculate 



°° ih 1 <- -» <- -> 



(42) 



with H = W(H) = \{p 2 + x 2 ). Since H is quadratic in both x and p, the sum 
in (|42|l only has to range from to 2. H * p = Ep therefore yields 



p 2 - ~r( d l + d D + iMxdp - P d x ) - 2E 



p = o 



(43) 



Separating (|43|l into its real and imaginary parts reveals two partial differential 
equations, 



(xd p -pd x )p = 



(44) 



and 



p 2 - (d 2 p +d 2 x )-2E 



P = o 



(45) 
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(|44l) shows that p(x,p) is a function only of x 2 + p 2 . It'll be more convenient 
to use u = 2(x 2 +p 2 )/h = 4H/h. We write p(x,p) = p(u), and substitute into 
(|45l) . The chain rule gives 

n2 d fdu d\ 4^ 16a;2 



and 



d 2 = \d u + l ^fd 2 u (47) 



similarly. Substituting into (|45|l gives 

= ^-d u -ud 2 u -^p{u) . (48) 

This is the differential equation to be solved to determine the Wigner function. 
It is no more difficult to solve than the Schrodinger equation in the operator 
formulation (but does not lead to Hermite polynomials, as we'll see). 

Setting p(u) = e~%L(u) simplifies considerations. Substituting into l|48|l 
gives 

= [udl + (1 - u)d u + L(u) . (49) 

We can look for a series solution, by substituting L{u) = J2^Lo ^a,u a . The 
resulting recursion relation is 

(a + l) 2 ^ a+1 = (a- f + ^4 ■ (50) 

For a normalizable solution, we need L(u) to be a polynomial. That is, t a must 
vanish for all a greater than some finite n. The recursion relation tells us this 
will happen if E = (n + ^ )U, for some non-negative integer n. These are exactly 
the SHO quantum energies (recall that we set u> = 1). 

The recursion relation yields solutions Lq = 1, L\ = 1 — u, L 2 = l — 2u + ^u 2 , 
etc. The normalization constant can be fixed by requiring J dxdpp — 1. The 
final general result is 

(—l) n 

Pn - '-e- 2 -^L n (AH/H) , (51) 

where L n denotes the nth Laguerre polynomial. 

The SHO can also be solved in algebraic fashion, as in operator quantum 
mechanics. Defining 

(x + ip) , a* = —= (x — ip) , (52) 



2h V2h 
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we can write H — h (a* * a + 5) . The ladder functions satisfy 

[a, a*]* = 1, [if, a*]* = h . (53) 

The SHO ground state is described by a Wigner function obeying 

a * po = pa* a* = . (54) 

The form of the zeroth Wigner function can be found directly from these equa- 
tions. The others are found using 

Pn = a* n *p *( a *y n . (55) 




Wigner 



p x 

Figure 1: SHO Wigner Function for n = 4 




Figures ^ and [21 depict the Wigner function for the n = 4 stationary state 
of the SHO. Because it only depends on u, it has the circular symmetry seen 
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in Figure ^ Figure |2| therefore gives the Wigner profile for any straight line in 
phase space passing through the origin. 

The classical SHO with this symmetry follows a circular trajectory in phase 
space, centered on the origin. The naive expectation might be then that the 
Wigner function, describing the corresponding quantum state, is a spread-out 
version of this; a single, circular ridge located above the corresponding classical 
phase path. The Figures show that only part of the Wigner function looks 
like that. There are oscillations in p, and it even goes into negative values. 
This means it is not a true probability distribution; it is instead called a quasi- 
probability distribution. 

This feature partly explains why deformation quantization is less popular 
than other formalisms of quantum mechanics. We stress again, however, that 
it reproduces all the predictions of the more familiar operator methods. For 
instance, Figure [3] plots the probability density V{x) against x (it could also be 
V(x) against p, because of symmetry), found using (|38(l . The curve is identical 
to that found from the state vector solving Schrodinger's equation. It is also true 
that expectation values of operators, such as X m P n , calculated using l]4Up. agree 
perfectly with those calculated in the more familiar formulations of quantum 
mechanics. 

6 Conclusion 

Deformation quantization is another way to do quantum mechanics. It has some 
strange features, such as a quasi-probability distribution, and an exotic way of 
multiplying functions, the ^-product. But it is perfectly consistent, and its pre- 
dictions agree with those made using other formulations of quantum mechanics. 
Furthermore, it is independent of the other formulations. On another planet, it 
might be the method of doing quantum mechanics discovered first! 11 

We believe that studying it, along with other quantum methods, deepens 
understanding of quantum physics. For example, writing FG might allow one 
to slip into thinking that the product of operators is similar to an ordinary 
product of functions. But the ^-product simulates the product of operators, in 
the way discussed above, and so eqn. makes clear that products of operators 
are tricky things. 

Deformation quantization is a bit of a well-kept secret, especially at the 
undergraduate level. We hope our brief discussion of it here can serve as a 
gentle introduction to the subject for upper-level undergraduates, and perhaps 
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others. 
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